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We  consider  N  particle  systems  whose  probability  distributions  obey  the 
master  equation.  For  these  systems,  we  derive  the  necessary  and  sufficient  conditions 
under  which  the  reduced  n  particle  (n  <  N)  probabilities  also  obey  master  equations 
and  under  which  the  Ursell  functions  decay  to  their  equ-librium  values  faster  than  the 
probability  distributions.  These  conditions  impose  restrictions  on  the  form  of  the 
transition  rate  matrix  and  thus  on  the  form  of  its  eigenfunctions.  We  first  consider 
systems  in  which  the  eigenfunctions  of  the  N  particle  transition  rate  matrix  are  com¬ 
pletely  factorized  and  demonstrate  that  for  such  systems  the  reduced  probabilities 
obey  master  equations  and  the  Ursell  functions  decay  rapidly  if  certain  additional 
conditions  are  imposed.  As  an  example  of  such  a  system  we  discuss  a  random  walk  of  N 
pairwise  interacting  walkers.  We  then  demonstrate  that  for  systems  whose  N  particle 
transition  matrix  can  be  written  as  a  sum  of  one  particle,  two  particle,  etc.  contri¬ 
butions,  and  for  which  the  reduced  probabilities  obey  master  equations,  the  reduced 
master  equtions  become  in  the  thermodynamic  limit  those  for  independent  particles,  whi< 
have  been  discussed  by  us  previously.  As  an  example  of  such  N  particle  systems  we 
discuss  the  relaxation  of  a  gas  of  interacting  harmonic  oscillators. 
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ABSTRACT 


Wi*  cc*!is i do r  N  particle  systems  whose  probability  distributions  obey 
the  master  equation.  For  these  systems,  we  derive  the  necessary  and  sufficient 
conditions  under  which  the  reduced  n  particle  (n  <  N)  probabilities  also  obey 
master  equations  and  under  which  the  Urscll  functions  decay  to  their  equili¬ 
brium  values  faster  than  the  probability  distributions.  These  conditions 
impose  restrictions  on  the  form  of  the  transition  rate  matrix  and  thus  on 
the  form  of  its  eigenfunctions.  We  first  consider  systems  in  which  the  eigen¬ 
functions  of  the  N  particle  transition  rate  matrix  are  completely  factorized 
and  demonstrate  that  for  such  systems  the  reduced  probabilities  obey  master 
equations  and  the  Ursell  functions  decay  rapidly  if  certain  additional  condi¬ 
tions  arc  imposed.  As  an  example  of  such  a  system  we  discuss  a^random  walk  of 
N  pairwise  interacting  walkers.  We  then  demonstrate  that  for  systems "whose  N 
particle  transition  matrix  can  be  written  as  a  sum  of  one  particle,  two 
particle,  etc.  contributions,  and  for  which  the  reduced  probabilities  obey 
master  equations,  the  reduced  master  equations  become  in  the  thermodynamic 
limit  those  for  independent  particles,  which  have  been  discussed  by  us  pre¬ 
viously.  As  an  example  of  such  N  particle  systems  we  discuss  the  relaxation  of 
a  gas  of  interacting  harmonic  oscillators. 


1 


I 


1.  lNTK'f.dCTiON 


In  our  previous  papers  in  this  series^ (hereafter  referred  to  as 
1,  11,  ill,  respectively)  on  the  decay  of  correlations  in  various  model  systems, 

we  hove  ceno latently  found  the  relations: 

Pn(u";t)  -  P*oJ(«*n)  *  PjCujt)  -  P}0j(c.);  n  a  1  (1.1) 

Un(an;t)  -  l^0)(un)  'vjl’1(a;t)  -  p{0)(a)]n  (1.2) 

where  Pn(an;t)  stands  for  the  n-particle  distribution  function  at  time  t  with 
dynamical  variables  a11  a  (a^  ,a2, . . .  ,an)  ;  Rn(an;t)  stands  for  the  n-particle 
relaxantsv  '  exemplified  by  the  Ursell  function,  Un,  in  papers  I  and  II  and 
the  C-function  in  paper  III;  P^  denotes  the  one  particle  distribution  function 
and  the  superscript  zero  denotes  the  equilibrium  value.  Equation  (1.1) 
describes  the  asymptotic  relaxation  of  the  n-particle  distribution  functions 
to  their  equilibrium  value  P ^  and  Eq.  (1.2)  describes  the  asymptotic  relaxation 
of  the  n-particle  relaxant  to  its  equilibrium  value  .  Equation  (1.2)  implies 
the  important  result  that  Pn(t),  n  >  1,  relaxes  to  a  functional  of  lower  order 
distribution  functions  £Pn_j(t),  Pn_2(t) » "  •  »pi  C^)j  as  |Pj(a;t)  -  p[°^(a)j  . 

These  remits  hold  for  such  diverse  systems  as  non-interacting,  initially 
correlated,  particles  in  contact  with  a  heat  bath  (I),  an  infinite  chain  of 
coupled  harmonic  oscillators  (II),  and  the  relaxation  of  spins  in  Glauber's^ 
one-dimensional  Ising  model  (III).  It  is  evident  that  there  must  be  an  under¬ 
lying  physical  basis  for  these  identical  results  on  the  decay  of  correlations. 

In  this  paper  we  study  the  general  conditions  under  which  the  n- 
particle  distribution  functions  and  the  Ursell  functions  will  be  of  the  form 
shown  in  Eqs.  (1.1)  and  (1.2).  We  limit  our  consideration  here  to  many  body 
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systems  whoso  time  development  is  governed  by  a  muster  equation  as  exemplified 
in  papers  1  and  111. 

In  Section  2  we  develop  eigenfunction  expansions  for  the  total  pro-* 

bability  and  for  the  reduced  probabilities  P  (t) ,  n  ■  1 , . . . ,N  -  1 .  We 

hs  *  — n 

ilr> t)>i«  f*nd  sufficient  conditions  on  the  form  of  the  eigenfunctions 

und  the  master  operator  under  which  the  reduced  probabilities  obey  master 
equations  and  under  which  the  Ursell  functions  relax  asymptotically  to  their 
equilibrium  values  faster  thau  the  reduced  probability  distributions  relax  to 
their  equilibrium  values. 

In  Section  3  we  consider  systems  for  which  the  eigenfunctions  of 
arc  completely  factorized  into  single  particle  functions.  We  show  that  for 
such  systems  the  reduced  probabilities  obey  master  equations  and  the  Ursell 
'  functions  decay  more  rapidly  than  the  reduced  probabilities. 

In  Section  4  wc  consider  a  model  of  N  interacting  random  walkers  as 
an  example  of  a  system  for  which  the  eigenfunctions  of  are  completely 
factorized.  For  the  particular  random  walk  chosen,  in  which  the  only  transi¬ 
tions  are  those  in  which  one  particle  takes  a  step  to  the  right  while  another 
particle  takes  a  step  to  the  left,  n  decays  asymptotically  at  the  same  rate 
as  Un>  In  the  thermodynamic  limit,  thp  reduced  master  equations  reduce  to 
those  for  independent  particles,  and  the  asymptotic  relaxation  of  and  P^ 
is  again  that  found  in  Paper  1,  and  given  by  Eqs.  (1.1)  and  (1.2). 

In  Section  5  we  generalize  the  results  of  Section  4  and  demonstrate 
that  for  systems  for  which  ^  can  be  written  as  a  sum  of  one  particle,  two 
particle,  etc.  contributions,  and  for  which  the  reduced  probabilities  obey 
master  equations,  the  reduced  master  equations  become  in  thermodynamic  limit 
those  for  independent  particles.  As  an  example  we  consider  in  some  detail 
a  system  of  N  interacting  harmonic  oscillator*  which  undergo  transitions  as 
result  of  two  body  resonant  collisions. 
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1. 1  GliNRJNCTI  ON  EXPANSIONS 


n 


Wc  consider  a  system  containing  N  particles  und  with  a  time  dependent 

N  N 

multivariate  probability  l’^(a  ;t) ,  where  u  ■  (aj  ,0^, . . . .  ,a^)  where*  is  some 
property  or  set  of  properties  of  the  i**'  particle.  We  assume  that  the  time 
dependence  of  P^,  is  governed  by  the  master  equation 

|t  PNC«Nit)  •  £  AN(aN,YN)  PN(YN;t)  (2.1) 

yN 

N  N 

where  the  transition  rate  A^(a  ,y  )  is  related  to  the  N  particle  "gain  and 
loss"  transition  rates  BN(aN,yN)  by 

^(a  ,y  )  «  BN(a  ,Y  )  -  4  N  N  2-  BN(p  ,a  )  (2.2). 

o  jjM 

Here  and  in  the  equations  to  follow  6  N  N  is  the  Kronecker  delta  for  discrete 

<*  »Y 

variables  and  the  Dirac  delta  function  for  continuum  variables.  We  shall  write 
the  solution  to  Eq.  (2.1)  in  terms  of  eigenfunction  expansions  and  investigate 
under  what  conditions  the  reduced  probabilities  Pn(an;t),  n  •  1,...,N  -  1  obey 
master  equations  of  the  form 


A  (a 
n 


n 


n. 
Y  ) 


PnCYn;t) 


(2.3) 


We  shall  also  investigate  the  conditions  under  which  the  Ursell  functions,^ 
Un(an;t),  decay  to  their  equilibrium  value  faster  than  the  Pn(an;t)  decay  to 
their  equilibrium  values. 

Equation  (2.1)  can  be  written  in  operator  notation  as 


St 


(2.4) 


with  the  formal  solution 
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At 

]'N(tJ  -  e"'  •  l'N(Oj 


(2.5) 


Tli*  rijlit  iiiul  lift  eigen function  ^  ^  and  K  of  obey  the  equations 


*  ~N , X  °  XN’  -^N.X 
’  M  N 


(2.6) 


(2.7) 


i»j  X  *  *Vl  X*  "  6X  A* 

n  ’*N 


(2.8) 


where  the  A^  arc  the  eigenvalues,  whose  spectrua  aay  be  discrete  and/or  con¬ 
tinuous.  We  consider  here  only  systeas  for  which  the  equiiibriua  distribution 
is  unique,  i.c.  «  0  is  a  nondegenerate  eigenvalue.  The  other  eigenvalues, 
however,  may  be  degenerate  and  the  eigenfunctions  should  have  a  label  denoting 
the  state  as  well  as  the  eigenvalue.  We  shall  suppress  this  label  for  ease 
of  notation,  bquation  (2.8)  is  understood  to  iaply  that  the  eigenfunctions 
for  different  states  aro  orthogonal.  The  probability  P^t)  can  now  be  written 

V«)  ■  E  •*"*  (±n,xn  •  *N(0))  VN  <j  »> 

where  the  symbol  ]£ denotes  the  sua  over  the  discrete  spectrua  and  the  integral 
over  the  continuous  spectrua  of  eigenvalues. 

From  detailed  balance 


where 


.  ,  N  N,  O(o),  N.  a  r  N  N.  D(o),  N. 
AjjCa  ,Y  )  P,}  '(Y  )  -  ^(y  ,o  )  P,J  "(a  ) 


P,S0)(aN)  m  PN(aN;-) 


(2.10) 


(2.11) 


is  the  equilibrium  probability  distribution,  and  froa  the  fact  that  is  a 
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stochastic  matrix,  it  fuiiow.s  that  the  eigenvalues  are  real  and  non  - 
pus i t  ive^  ami  that 

n0)‘“''')]’1  y»N‘«N>  (2-12)' 

If  *>  0  is  part  of  the  discrete  spectrum 


n(*) 


"  -^M,0 


(2.13). 


If  XN  >  0  is  part  of  the  continuous  spectrum 

P*o)  «  0  (2.14) 

and  in  liqs.  (2.10)  and  (2.12)  Q  should  bo  used  instead  of  . 

Tlic  reduced  probabilities  Pn(on;t),  n  ■  1,...,N  -  1,  for  the  n  particles, 
ij,..,in,  are  defined  by 


Pn(on;t:)  -  £  PN(aN;t)  (2.15) 

a"  n 

where  the  sum  is  over  the  variables  of  the  N-n  other  particles.  It  follows 
from  liq .  (2.9)  that 


•  £  «V(vn  •■V’)) 

No  now  investigate  the  conditions  under  which 
Eq.  (2.3),  i.e.  when  P  (t)  is  given  by 


(2.16) 


a 

Pn(an;t)  is  a  solution  to 


io(,> 


■  c 


At 

•Ml 


In<°> 


or 


%«•**) 


E 


V 


K  •  £"(0)) 


n,A 


(o  ) 


(2.17) 


(2.18) 
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Kiir  1 1  u  ,  aiul  I  t  arc  the  lijjlii  and  left  cl  (ten  functions,  and  A  the 
ii  '  n 

cii-i-nv.iliu-r.  of  A  .  It  is  obvious  that  the  A  must  form  o  subset  of  the  A„. 

•  II  N 


Clearly,  l.q.  (. 2 . 1 1» j  reduces  to  f:q.  12. 1R)  if  and  only  if,  for  each 


> ,  either 


ij  £  R..  .  uV-  R.  ,(«")"  0  (2.10) 

N-n  *'tAN  n*  N 

(I 

or 


(2.20). 


I.o t  us  consider  a  >N  for  which  Eq.  (2.19)  does  not  apply;  for  this  A^ 
liq.  (2.20)  must  apply  and  we  can  write 


u  . 

n.A 


(®  ) 


N 


(2.21) 


since  bq.  (2.20)  must  be  true  for  all  P^(0).  This  subset  of  A^  for  which 


hq.  (2.21)  applies  is  identical  to  the  set  of  Ar.  It  then  follows  from 
hqs .  (2.12),  (2.13)  and  (2.21)  that 


*,AN  R  (on) 

n,o 


R  i  <•") 
n,AN 


(2.22). 


The  physical  significance  of  the  fact  that  the  reduced  distribution 
functions  obey  the  mastor  equation,  Eq.  (2.3),  is  more  readily  understood 
from  the  condition  that  must  be  imposed  on  A^.  Tnis  condition,  which  is  of 
course  equivalent  to  Eq.  (2.21),  is 


An(on.vn)  »  £  An(oN,YN)  (2.23), 

o  ’n 

i.c.  the  transition  rates  for  tho  n  particle  subsystem  must  be  independent 
of  the  initial  states  of  the  other  N-n  particles.  Equation  (2.23)  imposes 
severe  restrictions  on  the  properties  of  A,,  but  these  conditions  are  not  severe 
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enough  lo  determine  uniquely  its  function:)!  lorm.  he  discuss  some  sufficient 
conditions  for  liq.  (2.23)  to  hold  in  Section  3. 

he  now  investigate  the  conditions  on  the  eigenfunctions  of  the  master 
operator  under  which  the  Urscll  functions  relax  to  their  equilibrium  form 
faster  than  the  probability  distributions  relax  to  their  equilibrium  form, 
he  start  with  the  discussion  of  £2  an<l  the))  extend  our  considerations  to 
LJ^.  The  two  particle  Urscll  function  is  defined  by 

U2(o1,m2;t)  «=  P2(a1#a2;t)  -  ^(a^t)  Pj(u2;t)  (2.24) 

with  the  property 


£  U2<Va2:t*  “  2  U2(al,a2:t)  "  0  (2.2S) . 

°1  “2 

The  time  dependence  of  U2  is  given  explicitly  by 

u2(Va2;t)  "  £  e  (^n,xn  R2,AN(or°2J 

N  (2.26) 

'  &  *  N  "  fa  '  V0))  fa  '  ■En<0))  "‘•>n(°i)  RW°2 


) 


where  we  have  used  Eq.  (2.16)  and 


U*  ) 


c 


E 

N-n 


a 

Equation  (2.26)  has  the  properties  of  Eq. 


(2.2S)  since 


(2.27). 


%  X  (°N>  ■  «X  0  l2-28) 

which  follows  from  Eqs.  (2.8)  and  (2.12).  If  no  restrictions  are  imposed  on 
the  eigenfunctions  FT.  the  asymptotic  decay  of  U.  is  the  same  as  the 
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asymploi  ie  demy  of  I',,  which  is  governed  by  A*  which  is  the  hi  sliest  non  zero 

eigenvalue  for  which  K.,  is  non  zero.  Tlio  asymptotic,  decay  of  P.  will  be 

the  Sana'  as  I  he  asymptotic  decay  of  P.,  if  II.  .  is  non  zero  for  the  highest 

*AN 

in.'ii  zero  e Lgcwaiue ;  if  R.  .  is  zero  for  this  eigenvalue,  P.  will  decay 

*  N 


In  order  for  to  decay  faster  than  is  necessary  that  the 

coefficient  of  the  term  zn  Eq.  (2.26)  containing  the  highest  non  zero  eigen¬ 
value  be  zero  for  all  initial  conditions.  This  will  be  the  case  if  and  only 
if 


**2  A°^al 
•  AN  1 


»a2)  = 


Rl,0(al)Rl.Xj(o2^ 


*  Rl,Ajtal)Rl,o'e,2) 


for  each  eigenstate  of  A°. 

The  general  expression  for  is  given  by 


(2.29) 


u„(cn;t) 


E  (-Dk(k-l)IPm  (a. 

£  1 


»a;  ;t) 


P  (a. 


m. 


m 


m-m^+1 


m 


-  E  (-Dk(k-l)l 

c 


(k) 

N 


A  (\x^ 


(2.30) 


where  the  sum  is  over  all  partitions  5  of  n  particles  in  subgroups,  k  is  the 
number  of  subgroups,  and  nn  the  number  of  particles  in  the  i^  subgroup.  The 
n  particle  Ursell  function  has  the  important  property 


E  u>n;t)  =  o 

o 1  n 


(2.31) 
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where  the  subscript  i  denotes  any  of  the  n  particles.  An  extension  of  the 
argument  of  the  last  paragraph  yields  the  result  that,  the  U^,  m=2,,..,n,  decay 
asyiapt ot i cally  faster  than  cxp(A°t)  if  und  only  if 


(2.32). 


Note  that  in  this  case  the  time  dependence  of  _Un  is  independent  of  whether  or 
not  the  equilibrium  probability  distribution^0^  factorizes  into  single 
particle  functions. 

If  £n(t),  n=l,...,N-l,  is  a  solution  of  the  reduced  master  equation, 
Eq.  (2.3),  so  that  Eqs .  (2.19)  through  (2.22)  apply,  the  Urscll  functions 
will  not  decay  any  faster  than  the  probabilities  unless  additional  assumptions 
are  made.  Sufficient  additional  assumptions  arc: 

i)  The  equilibrium  distribution  is  factorized  into  single 
particle  functions,  i.e. 


P*0) («")  -  17  Pj00^.)  ; 


(2.33) 


ii)  The  eigenfunction  R.  '*(«.)  /  0  for  each  of  the  eigen* 

1,AN  1 

states  of  for  one  of  the  particles  of  the  set  n. 

Under  these  conditions  it  follows  from  Eq.  (2.22)  that  for  each  eigenstate 

of 


R  ,(o.)  *  0  ,  j  i  i  (2.34) 

1  ’  N  J 

and  thus  Eq.  (2.22)  is  identical  with  Eq.  (2.32).  Therefore  U^,  m=2,...,n, 
decays  asymptotically  faster  than  P  ,  m=l,...,n. 
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3.  COM!*! !  n  i.Y  FACfORIZEL'  lilULNl-lWCTIONS 


As  we  h;« v f  shown  in  the  last  section,  the  specification  of  the 
necessary  conditions  for  the  fast  decay  of  the  Ursell  functions  is  quite 
complicated.  In  this  section  we  shall  consider  some  fairly  stringent 
.sufficient  conditions  for  tno  fast  decay  of  the  Ursell  functions  tnd  for 
the  validity  of  the  master  equation  for  the  reduced  distribution  functions, 
hq.  (2.3). 

We  assume  that  all  of  the  eigenfunctions  of  the  master  operator  AN 
arc  completely  factorized  into  single  particle  functions,  i.e. 

N  N 

R  (a  )  =  I!  R  (a.)  (3.1) 

N,AN  i=l  Vi  1 

where  wc  use  the  notation 


N, 


=  ^N^V1 *  * ' ’ * VN^  *  AN^V  ^ 


(3.2) 


with  XN(0,...,0)  =  0.  Under  these  conditions  Eq.  (2.9)  can  be  written 


where 


N  v  Vw)»/  \  N 

V“  i*)  ■  £  «  •  4(0)j  .qx  \<ai> 

N  N  /  N 

L..  1  (a  )  =  n  R  (a. )  /  R  (a. )  s  [J  L  (o'.) 

N*V  i-l  V  1  /  0  1  i-1  V  1 


(3.3) 


(3.4) 


.  n(0) 


and  R0(au)  =  (cu)  is  the  one  particle  equilibrium  distribution.  Equation 
(3.4)  follows  from  Eq.  (2.12).  The  orthonormality  condition,  Eq.  (2.8),  is 
equivalent  to 


x'  x 


L  Lv.(ai)  Rv!l°i)  *  4v.,v! 
1  1  ’  ’  ’ 

If  we  take  =  0  this  yields 


(3.5)  . 
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(3.6) 


ER  ,(o.)  =  6  , 

v !  j  o,v! 


u.  l 
l 


liquation  (2.16)  for  the  reduced  probabilities  now  becomes 

Xw(vN)t  /  \  n  N 

!>  Cr  n.«--\  -  V*  />  N  f  \  .  i>  rnti  rr  r>  ■)  77  f 

-  2*  (V>N  hr-  £  v'-V  pml  ».v, 

(3.7) 

where  we  have  used  liq.  (3.6).  The  sum  over  Vn+]>***»VN  Ecl-  (3.7)  can  be 
immediately  performed  to  yield 


A  (vn)t  /  \  n 

Pn(*n;t)  -  2  en  L  •  4(0)1  .17  Rv.(-) 

n  \  n  /  i=l  i 


(3.8) 


where 


*n(vn)  s  AN(vn,o,.. .,o) 


(3.9) 


and 


L  (an)  «  /7  L  (a.)  (3.10). 

’  n  i=l  i 

liquation  (3.8)  is  obviously  of  the  form  of  Eq.  (2.18)  and  thus  the  reduced 

probability^  obeys  the  master  equation,  Eq.  (2.3).  Therefore  the  factori- 

zation  of  the  eigenfunctions  FT,  .  as  given  in  Eq.  (3.1)  is  a  sufficient 

- "^,AN  - 

condition  for  the  validity  of  Eq.  (2.3). 

We  shall  now  investigate  the  properties  of  the  Ursell  functions  when 

all  the  eigenfunctions  factorize.  It  follows  from  (3.8)  and  (3.10)  that  the 

terms  in  Eq.  (3.8)  for  which  =  0,  j  =  l,...,n,  are  of  the  form 

Pn_j(an’*;t)  p|°^(Oj).  Thus  Eq.  (3.8)  can  be  rewritten 

n 

(3.11) 
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where  the  sum  over  n  is  over  all  partitions  of  n  particles  into  two  sub¬ 
groups  where  m  goes  from  1  to  n.  The  notation  '  implies  the  sum  over 

v* 

all  possible  non-zero  values  of  Substitution  of  Eq.  (3.11) 

into  Eq.  (2.30)  for  U^(t)  yields 

m  m.  \ 

Am  (v  1 )  A  (v  K)  It 

ml  \  / 

(3.12) 


Un(an;t) 


■  £  (-u1 
£ 


(k-l)l 


£ 

n 


exp 


( 


x 


W°j 


n 

n 

i»l 


where  we  have  used  the  fact  that  factorized  probabilities  do  not  contribute 

to  the  Ursell  functions.  Note  that  because  of  the  restriction  on  the  sum 

over  vn.  all  the  A  that  appears  in  Eq.  (3.12)  are  non  zero.  The  asymptotic 
"i  m.  m^ 

decay  of  is  determined  by  the  maximum  value  of  Am  (y  A)  A^^v 

vr...,vn  +  0. 

The  stringent  conditions  on  the  eigenfunctions  in  Sections  2  and  3  for 
the  rapid  decay  of  the  Ursell  function  LL(t)  do  not  provide  much  physical  insight. 
It  would  be  preferable  if  these  conditions  could  be  stated  in  terms  of 
properties  of  the  operators  and  ^  of  the  N  particle  and  reduced  n  particle 
master  equations.  We  have  been  unable  to  do  that  in  general.  In  Sections  4 
and  5  we  present  and  discuss  some  sufficient  conditions  on  the  and  ^  for 
the  rapid  decay  of  correlations.  These  conditions  do  yield  physical  insight 
into  the  types  of  interactions  which  lead  to  a  rapid  decay  of  correlations 
in  stochastic  processes  governed  by  a  master  equation. 
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4.  RANDOM  WALK  WITH  INTERACTIONS 


In  this  section  we  consider  N  pairwise  interacting  random  walkers  on 
a  finite  one-diinensional  lattice  with  periodic  boundary  conditions.  The 
sites  on  the  lattice  might  represent,  for  example,  the  internal  states  of  gas 
particles.  The  steps  of  the  random  walkers  on  the  lattice  mirrors  the  re¬ 
sults  of  the  binary  collisions  of  the  gas  particles.  The  features  of  the 
problem  which  are  of  particular  interest  to  us  are  the  asymptotic  decay  of 
the  Ursell  functions  and  the  structure  of  the  transition  rate  matrix  of  the 
master  equation  in  the  thermodynamic  limit. 

We  consider  N  random  walkers  with  binary  interactions  on  a  finite 
one-dimensional  lattice  with  l  sites  and  periodic  boundary  conditions.  The 
position  of  particle  i  on  the  lattice  is  denoted  by  nK»l,...,fc;  with 
i=l,...,N.  In  our  model,  any  site  m  can  be  occupied  by  more  than  one  particle 
subject  to  the  restriction  that  the  total  number  of  particles  equals  N.  We 
write  the  master  equation  of  the  N  particle  system  as 


d 

5t 


N 


(mN;t) 


,m^-l , . . . ,m^+l , . . . ,m^;t)  -  P^(m 

(4.1). 


where  a  is  a  constant  transition  rate  for  the  transitions  jik  — ►  m^-1; 

dk  — ►  nij+1 .  This  equation  describes  a  system  with  two  particle  interactions 

such  that  one  particle  uioves  one  step  to  the  right  while  the  other  particle 

moves  one  step  to  the  left  while  all  other  particles  remain  in  the  same 

N 

state.  In  this  model  Y*  m.  is  conserved  in  the  sense  that  for  all  states 

fai  1 

that  occur  in  Eq.  (4.1), 


C 


♦ 


kt 


(4.2) 


where  C  is  a  constant  and  k  is  an  integer. 
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The  eigenfunctions  of  the  master  operator  are 


_  ,  N.  a-N  f  2ui  N  Nl 

N  ym)  exp  v  m  j 

N 

=  n  R,  (m.)  ,  v.  «  0, 

j'i  vj  ’  3 

where 


(4.3) 


R  (m.)  =  t"1  exp  v  m  (4.4). 

vj  J  1  J  J 

The  eigenvalues  are 

V“N)  ■  ir  2  (vi  -  V  - ']  (4-5)- 

1 

It  is  clear  that 
yvN)  -  o 

if  and  only  if 

vk  »  v,  k=l,...,N  (4.6) 

for  all  vk=0, . . . ,4-1.  The  zero  eigenvalue  &N  is  l  fold  degenerate  which  is 
a  consequence  of  the  conservation  law  expressed  in  Eq.  (4.2).  The  highest 
non-zero  eigenvalue  for  and  is 

*  XN^V  +  v*  v'-***v)  ’  jr  (N-l)(cos  j—  -  1)  *  -  A  (4.7). 


The  N  particle  eigenfunctions  in  this  model  are  completely  factorized 
into  ;ne  particle  functions  and  thus  we  are  dealing  with  a  special  case  of 
the  systems  discussed  in  Section  3.  The  reduced  probabilities  obey  master 
equations  of  the  form  of  Eq.  (2.3).  The  reduced  eigenfunctions  are  given  by 
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(4.8) 


n 

R  (mn)  -  n  Rv  («J 

n,v  j«l  j  1 


and  the  reduced  eigenvalues  are 


An(v")  =  XN(vn,0,...,0) 


2o(N-n) 

N 


S  [“tv1]  |  [«r<vV-‘] 

1 


(4.9). 


It  should  be  noted  that 


Xn(v")  =  0;  n=l , . . . ,N-1 


(4.10) 


if  and  only  if 


v^«0,  k  »  1,  2, . . . ,n 


(4.11) 


Thus  the  zero  eigenvalue  for  all  reduced  probabilities  is  non  degenerate. 
The  highest  non- zero  eigenvalue  X*  for  Pn(mn)  is 


n 


xn  ■  "  r  r  - 1}  ■  • A 


(4.12) 


and  thus  all  the  reduced  probabilities  decay  asymptotically  to  their  equili¬ 
brium  values  as 


Pn(mn;t)  -  P<o)(mn)  %  cne’At 


(4.13) 


where  cn  is  a  constant  which  depends  on  the  initial  conditions  and  where 


p£V)  -  *-n 


(4.14) 


The  expression  for  Un(mn;t)  is  given  in  Eq.  (3.12).  The  highest  non¬ 
zero  eigenvalue  appearing  in  this  expression  is 
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(4.15) 


2an(N-n) 

N 


/  2n  ..  N-n  . 
(cos  j —  1)  ■  -  jjTy  nA 


for  n>l,...,N-l.  Thus,  Un(mn;t),  n«2,...,N-l,  decays  asymptotically  to 
its  equilibrium  value  *  0  as 

.  E^IllAt 

Un(mn;t)  ^  e  *  (4.16) 


where  c^  is  a  constant  which  depends  on  the  initial  condition.  Note  that 
Uj (m; t)  h  Pj (m;t) . 

From  the  form  of  A*  in  Eq.  (4.7)  it  follows  that  the  N-particle 
functions  ^(t)  and  U^(t)  decay  asymptotically  to  their  equilibrium  forms  as 

4(0  -  P*0’  '  cNe'At  (4.17) 

and 

^(t)  -  U<0)  '  Cj*e‘At  (4.18) 

where 

^(“0  -  40)  -  i"N  (4.19) 

and  depends  on  the  initial  conditions  because  of  the  degenerate  zero 
eigenvalue  as  given  in  Eq.  (4.6). 

We  now  discuss  the  behavior  of  the  Ursell  functions  U  (t)  for  finite 

~n 

N  and  in  the  thermodynamic  limit  as  N— *•.  For  N  finite,  n(t)  decays 

asymptotically  at  the  same  rate  as  U  (t),  which  follows  immediately  from 

Eq.  (4.16).  The  fastest  decaying  Ursell  functions  are  those  for  which 
N 

n  dr  2*.  The  slowest  decaying  Ursell  functions  are  those  with  n«N  and  n«N. 
The  Ursell  functions  for  nacN  decay  slowly  owing  to  the  conservation  law 
given  in  Eq.  (4.2).  In  the  thermodynamic  limit,  the  asymptotic  decay  of  the 
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reduced  probabilities  P  (t)  is  still  given  by  Eq.  (4.13)  while  tho  decay 

— n 

of  the  Urscll  functions  U  (t)  is  given  by 


v*> 


c* 

n 


-nAt 


(4.20) 


which  follows  from  liq.  (4.16).  In  the  thermodynamic  limit,  our  model  systems 
thus  has  the  identical  asymptotic  behavior  for  £n(t)  and  U^t)  as  previously 
found  in  1  and  II. 

The  following  question  immediately  arrises:  why  does  this  model 
system  with  interactions  behave,  in  the  thermodynamic  limit,  in  the  same  way 
as  the  independent  particle  systems  considered  in  1?  The  operator  in  the 
master  equation  (4.1)  can  be  written  as 


V"  *1  >  ■  !  |  j[sVMi  VMi  MM  4Vk]- 


(4.21). 


The  operator  for  the  reduced  master  equation  is 
A. 


I  .  %  t  I  «  ,  n  I  1-  t  ) 

n  n  m  jl  mm  v«k| 

1  1 

S  |[(S"k'1-<'i  *  V*'0  k?i  2V,q"| 


(4.22) . 

In  the  thermodynamic  limit,  ^  becomes 

Vqk 

(4.23) . 


An(»n.qn) 


£  [V 


Mi 


♦  6 


v1* 


26. 


Vqi 


] 


n 

n 

Mi 

1 


.  18  . 


Thus,  in  the  thermodynamic  limit,  the  n-particlc  transition 
sum  of  single  particle  transition  rates  and  is  identical  in 
of  Paper  I  for  independent  particle  dynamics. 


rate  A  is  the 
form  to  Eq.  (2.15) 


-  19  - 


5.  MASTER  EQUATION  MODELS  FOR  INTERACTING  GASES 


In  this  section  we  consider  systems  containing  N  identical  interacting 
particles  in  which  the  operator  A^  for  an  N  particle  master  equation  can  be 
written  as  a  sum  of  one  particle,  two  particle,  etc.  contributions.  We  show 
that  in  the  thcrmoJynamic  limit  as  N— ♦»,  the  reduced  master  equation  operators 
An  can  be  written  as  a  sum  of  one  particle  terms  as  long  as  Hq.  (2.23)  applies. 
We  illustrate  our  discussion  by  treating  a  system  of  N  interacting  harmonic 
oscillators . 

We  consider  an  N-particle  transition  rate  matrix  of  the  form 


.  ,  N  N\ 
An(o  ,y  ) 


N  N 

E  a  n 

i=l  1  11  k/i 

1 


6  + 
VYk 


]_ 

2N 


(a.  ,a  JY .  ,Y  .  ) 


N 

n 

Wi  J 

1 


6  + 

VYk 


(5.1) 


5? 


N  (3)  N 


A  (tt.,a.,a  ;Y.,if.,Y.)  fl  6„  v  +  ••• 

6N‘  3  1  J  1  1  3  1  kj<i ,  j  ,1  VYk 

1  1 


where  A^ ,  Aj3^ ,  etc.  are  symmetric  when  the  particle  indices  are  permuted. 
Here  a|^  is  a  one  particle  transition  rate  matrix  which  has  the  property 


E  *i(1)!»iiV0  • 


(5.2) 


a. 

l 


(2) 

A£  is  3  two  particle  transition  rate  matrix  with  the  property 


E  a22,(vwV  3  0  ; 


(5.3) 


a.  ,a . 

i  J 


Aj3^  is  a  three  particle  transition  rate  matrix  with  the  property 
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etc.  The  form  of  in  Eq.  (5.1)  implies  that  the  systems  under  consideration 
can  be  described  in  terms  of  one  particle,  two  particle,  three  particle,  etc. 
interactions  as  is  the  case,  for  instance,  for  weakly  interacting  systems 
with  short  range  forces.  The  factors  N"^"^  in  front  of  the  sums  take 
account  of  the  frequencies  of  i-particlc  collisions.  Thus,  the  are  inde¬ 
pendent  of  N.  We  assume  that  for  the  operator  given  in  Eq.  (5.1),  the 

condition  (2.25)  holds  in  the  thermodynamic  limit  as  N— ►».  In  this  case, 

A  will  be  given  by 


.  ,  n  n, 
An(o  ;y  ) 


£  VvTi> 


n 


n 

Mi 

i 


6 

VYk 


(5.5) 


where 


A  =  Ad)  +  a(2)  +  a(3)  + 


(5.6) 


and  where 


A^1^(oi;yi)  =  A^^a.jy.) 


,(3) 


(5.7) 


(5.8) 


(»j,Y.)  ■  lim  ±  E  E  a'2)(c.  ;y  ) 

N— j  =  l  a.  J  J 

m  3 

(c  tp-lim  4  E  E  Z  A33)(oi,“j'“»:yi,Yj'V  (S-9) 

N— mo  N  j  ,£=1  a .  a  J  J 


etc.  In  writing  Eqs.  (5.8)  and  (5.9)  we  have  assumed  that  the  two  particle, 
three  particle,  etc.  contributions  to  A^  obey  the  condition  (2.23)  in  the 

thermodynamic  limit.  In  many  systems  A^  (a.  ,a .  ;y.  ,y.)  will  be  independent 

a .  J  1  J 

J 


of  Yjj  in  other  systems,  as  we  shall  see  below,  it  is  necessary  to  take  the 
sum  over  j  to  obtain  a  quantity  independent  of  in  thermodynamic  limit. 

Note  that  if  An  has  the  form  given  in  Eq.  (5.5),  the  results  given 
in  1  for  independent  particle  systems  will  again  apply. 

We  now  consider  the  example  of  a  gas  of  interacting  harmonic  oscil¬ 
lators  which  exchange  only  vibrational  energy.  The  N  particle  master  equation 

(7) 

has  the  form  ' 

d  n  ,  N  V  .  ,  N  N,  _  ,  N  ^ 

PN(m  ;t)  =  L,  A^(m  ,q  )  PN(q  ; t)  = 

N 

q 


^ ]  ^>1  (m j  + 1 )  P^  (m^ , . . . , nr  - 1 , . . .  ,m j  + 1 , .  •  •  ,m^  j  t) 

1 


(5.10) 


(®i+1  PN(mN;t)j 


where  £  is  a  rate  coefficient  and  where  jik=0,  1,...  with  i=l,...,N,  denotes, 
the  states  of  the  ith  oscillator  with  energy  E(nh)  =  nrhui.  The  dimensionless 
mean  vibrational  energy  per  particle,  e  is  given  by 


N 


e  ■  Nib  2  B(»P  *  jj  Z 


m. 


(5.11). 


m. 

l 


m.  =1 
x 


Clearly,  this  mean  energy  is  a  constant  which  is  conserved  during  the  relaxa¬ 
tion  process.  In  order  for  the  relaxation  process  to  take  place  on  the 

constant  energy  shell  with  energy  Ne  it  is  necessary  to  choose  initial  con- 

N 

ditions  such  that  PN(m  ;0)  is  non-zero  only  when  ^  m.  =  Nc.  The  N-particle 

in. 

master  operator  has  the  form  1 

ANCmN;qN)  *  ^  £  A^2)(m  ,m  ;q  q  )  /7  6 

i/j  1  J  1  J  k=1  Vqk 

1  Mi,j 

where 


(5.12) 


22  - 


A,2^  (m.  ,m.  ;q.  ,q.)  =  a  <  m.  (m.  +  l)  6  .  6  . 

2  i*  l  i  J  n'i-l,qi  nu+l,q.. 

CS.13). 

♦  (m.+l)  in.  6  ,  5  ,  -  f  (m .  + 1 )  m.  +  m. (m.  +  l)  ]  6  6  } 

i  J  m.+l.q.  l  }  l  3  J  Vqi  ^.q.  f 

liquation  (5.12)  is  clearly  a  special  case  of  Eq.  (5.1).  The  two  particle 

(2) 

transition  rate  matrix  has  the  properties 

C  (mi  »ll j )  5  J]  A^2)(mi,mj;q.,qj)  = 


m. 

J 


aq .  fin.  6  .  ♦  (m.+l)  5  .  -  (2m. +1)  6  l  (5.14) 

l  i  mi“l*c*i  1  V1,qi  1  nii  ,qi  / 


+  a 


f (m. +1)  6  .  -  m.  6  l 

i  i  "‘i*1  »q±  i  vv 


and 


A(2)  (m.;q.) 


1  & 

lim  K  £  C(m.;q  q  ) 
N-*»  j*l  1  1  J 


(5.15) 


m.  6  , 

l  m.-l,q. 
i  ’  u 


+  (e+1)  (m.+l) 


mi+i *qi 


(2emi 


+  e  + 


mi) 


m 


where  the  limit  N— +*  is  taken  keeping  e  fixed.  Equation  (5.15)  leads  to  an 
(2) 

Av  *  of  the  form  given  in  Eq.  (5.8).  From  the  relation  given  in  Eq.  (5.6) 
it  then  follows  that  A^  is  of  the  form  of  Eq.  (5.5)  with  the  n  particle  master 
operator  given  as  a  sum  of  1  particle  transition  rates.  The  results  of  Paper  I 
are  thus  directly  applicable  to  the  system  of  interacting  oscillators.  The 
eigenfunctions  and  eigenvalues  of  this  one  particle  transition  rate  matrix 
Aj Cai »Yi J  can  be  found  in  the  paper  by  Montroll  and  Shuler. ^ 

The  results  of  this  section  can  be  briefly  summarized  as  follows.  If 
the  transition  rate  matrix  A^.  for  the  N-particle  master  equation  can  be  written 
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in  the  form  of  Eq.  (5.1),  which  is  the  case,  for  instance,  for  systems  with 
weak  interactions  and  short  range  forces,  and  if  the  operator  is  of  the 

form  given  in  Eq.  (2.23),  i.e.  if  the  N-particle  master  equation  reduces 
to  an  n-particlc  master  equation  by  summation  over  the  other  N-n  particles, 
then  in  the  thermodynamic  limit,  reduces  to  a  sum  of  1-particle  transition 
rates  A^ .  This  implies  that  for  all  such  systems 

V0  ■  ^0)  *  pi(t)  ■  £i(o)  ]n  (5.19) 

£n(t)  -  j^°3  *  ^(t)  -  Pf03]  (5.20). 

These  results  are  identical  with  those  obtained  by  us  in  1  for  systems  des¬ 
cribed  by  independent  particle  dynamics  which  are  also  characterized  by  1- 
particle  transition  rates  A  . 
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